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ABSTRACT 



^v^j | Context. Recentely, an increasing number of scientific publications making use of images obtained with near-infrared long-baseline 

interferometry have been produced. The technique has reached, at last, a technical maturity level that opens new avenues for numer- 
ous astrophysical topics requiring milli-arc-second model-independent imaging. The Very Large Telescope Interferometer (VLTI) 
will soon be equipped with instruments able to combine between four and six telescopes. 

Aims. In the framework of the VLTI second generation instruments Gravity and VSI, we propose a new beam combining concept 
using integrated optics (IO) technologies with a novel ABCD-like fringe encoding scheme. Our goal is to demonstrate that IO-based 
. combinations bring considerable advantages in terms of instrumental design and performance. We therefore aim at giving a full char- 

acterization of an IO beam combiner in order to establish its performance and check its compliance with the specifications of an 
imaging instrument. 

Methods. For this purpose, prototype IO beam combiners have been manufactured and laboratory measurements were made in the 
H band with a dedicated testbed, simulating a four-telescope interferometer. We studied the beam combiners through the analysis of 
throughput, instrumental visibilities, phases and closure phases in wide band as well as with spectral dispersion. Study of the polar- 
ization properties was also carried out. 

Results. We obtain competitive throughput (65%), high and stable instrumental contrasts (from 80% in wide band up to 100%+ 1% 
with spectral dispersion), stable but non-zero closure phases (e.g. 115° ±2°) which we attribute to internal optical path differences 
(OPD) that can be calibrated. We validate a new static and an achromatic phase shifting IO function close to the nominal 90° value 
£Sj ' (e.g. 80°±1°). All these observables show limited chromaticity over the H band range. 

■^J" , Conclusions. Our results demonstrate that such ABCD-like beam combiners are particularly well suited for interferometric combi- 

. nation of multiple beams to achieve aperture synthesis imaging. This opens the way to extending this technique to all near infrared 

■ wavelengths and in particular, the K band. 

£NJ ■ Key words, optical interferometry - integrated optics 
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^ resolution. The study of stellar surfaces, evolved stars, young 

. - - i stars, our galactic center and the heart of active galactic nu- Th e Very Large Tele s cope Interferometer (VLTI, (Scholler 

clei require access to direct imaging. Until now, a large frac- 2007; H aguenauer et alj 120081) ) will be equipped in 2012- 

tion of observations in the near infrared (NIR) were obtained 2015 with two seco nd-generati on instruments: Gravity 

with 2 to 3-telescope arrays, with little spatial frequency cov- (lEisenhauer et alJ |2008) and VSI dMalbet et alJ 120061 [2008) 

erage (so-called uv coverage), restricting the astrophysical in- which will be capable of exploiting the imaging capability of 

terpretation to a parametric one in most of the cases. However, the array by combining four beams for the first and six for the 

discriminating between different successful scenarios of com- second. The stringent requirements for these two instruments 

plex or rapidly-changing objects raises the need for images have triggered the interest in using integrated optics (IO) 

as model-independent as possible. This translates into the re- as a core technology for the beam combining function. The 

quirement to use as many telescopes as possible in order to ability to integrate a singlemode circuit on a substrate, able 

fill the uv plane and allow an unambiguous image reconstruc- to interfere all the beams offers numerous advantages both in 

tion. Until very recently, most of the images produced with terms of performance and ease of operation. Single-mode beam 

optical long baseline interferometers had moderate complex- combiners provide natural modal filtering, which associated 

ity and therefore did not bring additional information with re- with proper photometric calibration has been shown to lead to 

spect to parametric modelling. In our opinion, the complex- accurate visibility measurements. The compactness of the chip 

ity barrier where the reconstructed image adds meaningful sci- allows the instrument footprint to be minimized and the thermal 

entific value to the astrophysical interpretation was recently control to be optimized (further enhancing the calibration 

accuracy). No alignment is required, other than the injection 

Send offprint requests to: Myriam.Benisty@obs.ujf-grenoble.fr in the input guides, even though the combination scheme is 



1 . Introduction passed by lMonnier et all d2007l) and lZhao etaf] d2008l) using the 

MIRC instrument, an image plane 4-beam combiner using sin- 
Optical long baseline interferometry offers a unique way to di- g i e mode fibers , at the CHA RA interferometer (MonnierSl 
rectly probe astrophysical environments with milli-arcsecond 2006b; ten Brumm elaar et alll2l)05h . 
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Fig. 1. Upper panel: theoretical design of the integrated optics 4-way beam combiner allowing pairwise combination and using 
phase-shifting devices to produce 4 outputs in quadrature. We refer to each output using the index m,l,k. ml k is the k'' 1 output out of 
4, resulting from the combination of the beam m and I. The lower panel is a picture of a prototype that is 80mm long and 8mm wide. 



complex. Finally, this technology offers the flexibility to easily 
switch beam combiners to adapt to a particular situation (e.g. 
target, number of telescopes). 

Since the initial proposition by iKern et all dl996l) . LAOG 
and its industrial partner LETI/CEA have been developing the 
use of IO technology to interferometrically combine light beams 
in optic al waveguides lying on a solid substrate of a few cen - 
timeter (IKern et dJl99ftlMalbet et al.ll999l:lBe"rger et al.l2000l) . 
This instrumental research program has consisted in designing, 
fabricating and characterizing all the IO building blocks required 
to build an astronomical interferometric beam combiner. Several 
beam combining schemes have been implemented and tested. 
Some of them have led to successful on-sky demonstrations such 
as the VINCI/VLTI (2 telescopes) and IONIC 3/IOTA (3 tele- 
scopes) instrumen ts (Bergeret al. 2003; Le Bouquin et al.ll2004t 
Krausetal .1 120051: iMonnier et alJl2006ah . 

In the context o f VLTI second-generation instrument studies, 
iLeBouquinl (120051) have studied the global efficiency of a great 
variety of IO beam combiners. This study has concluded that one 
of the most efficient ways to combine four beams (e.g. 4 UT or 4 
AT) was to use a so-called "pairwi se static ABCD" schem e (in- 
spired by the visibility estimator of Shao & Staelin] d!977l) ). This 
IO circuit allows one to extract simultaneously four phase states 
of the coherent signal ind ependently for each of the six base- 
lines. We fabricated them dLabevdfeOOSl) . and in this paper, we 
present these new 4-beam combiners together with their com- 
plete laboratory characterization. They are probably the most 
sophisticated astronomical beam combiners built to date. The 
paper is organized as follow: in Sect. 2, the technology and the 
specific design of the beam combiners are described. In Sect. 3, 
we present the laboratory set up as well as the experimental pro- 
cedure; the characterization results are given in Sect. 4 and dis- 
cussed in Sect. 5. 



2. The beam combiner: technology and design 

Prior to fabrication, the IO circuit was designed and numeri- 
cal computation simulating the propagation of an electromag- 
netic signal was carried out to determine the expected proper- 
ties in terms of flux routing. Each IO function was checked and 
its throughput and flux distribution were optimized numerically. 



Y junction Phase shifting device 




(P« 2 = <pn' * 7T 
<p12 3 = <P12' + TT/2 

<ft2 i = qH^* n = qii2 2 + ir/2 



Fig. 2. Details of the beam combining function: for each inter- 
ferometric pair (e.g. [12]), one arm is shifted by 90° leading to 
four outputs in quadrature (with phases written as <pL to ^ 2 ). 
Combinations of beams occur in couplers that present two out- 
puts in phase opposition to maintain energy conservation. By 
recording the four phase states (ABCD-like, see the right fig- 
ure), one can retrieve the interferometric observables (amplitude 
and phase of the fringes). 

This step done, the simulation parameters were turned into tech- 
nological parameters and a photolithographic mask was fabri- 
cated. 

LETI uses a silica-on-silicon technology to fabricate IO cir- 
cuits. This technological process requires several photolitho- 
graphic steps to etch different layers. The beam combiners are 
made by depositing alternatively 3 doped silica layers on a 
silicon substrate. The second layer is etched to define chan- 
nel waveguides and the other two layers constitute the optical 
cladding. For the first time, the etching technolog y allows us to 
comp letely isolate each waveguide from the others ( Labev e"et al.l 
2006). The produced beam combiners have been designed to op- 
erate in the atmospheric H band and more recently in the K band. 

The so-called "pairwise static ABCD" beam combiner can 
be described as follows. Each beam combiner is designed to 
have 4 inputs and 24 outputs, allowing 6 interferometric pair- 
wise combinations, each one producing 4 phase-shifted outputs 
with a phase difference of 90°. For each injected beam, the light 
propagates through waveguides and is split in three in a tricou- 
pler (item (a) in Fig. [TJ to enter the combining function (consti- 
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tuted of Y-junctions and couplers). The light is then divided in 
two in a Y junction that acts like a classical beamsplitter (item 
b), each beam later being combined in a coupler (item c) with a 
beam coming from another telescope. A coupler allows a con- 
trolled power transfer between one waveguide and another. As a 
consequence of energy conservation, each coupler has two out- 
puts in phase opposition. In only one of these four arms, there 
is a phase-shifting device designed to change the phase of the 
propagating beam by 90° (Fig. [2]). This leads to four output 
beams, two of them being in phase opposition with an additional 
phase-shift of 90° with respect to the other two (e.g. tpyi and 
ipvi+ir, ipn+7r/2 and (pi2+7T+n/2, following Fig. [2]). The phase- 
shifting function is based on the variation of the effective index 
(i.e. index seen by the fundamental mod e propagatin g into the 
waveguide) with the waveguide diameter dLabevel2008l) . To cre- 
ate a phase shift, enlarging one of the two waveguides creates 
a difference in the effective index and leads to an optical path 
difference between two parallel waveguides of the same phys- 
ical length. In order to achieve an achromatic phase shift, the 
wavelength dependence is compensated by concatenating a few 
waveguide segments of different diameters separated by tapers 
(i.e. adiabatic functions) to avoid any loss due to discontinuities 
(Fig. |3). Since the photometry is extracted from a linear com- 
bination of the interferometric signal itself, the beam combiners 
have no dedicated photometric channels. This allows us to effi- 
ciently use all photons for the interferometric combinations. By 
design, each interferometric pair simultaneously gives access to 
four phase states in quadrature (ABCD-like but without tempo- 
ral modulation). These 4 measurements allow the visibility am- 
plitude an d phase to be r etrieved using the ABCD method de- 
scribed in iColavital d!999f) . In practice, the departure from ideal 
quadrature forbids the use of simple algorithms and leads us to 
consider a generalized algorithm capable of handling a realistic 
description of the beam combiner properties. 

Throughout the paper, the outputs are identified with the in- 
dex m, I, k, such as ml k , where m, I are the interfering beams, and 
fc=[1..4], the output for this combination (similarly, the A-C-B- 
D measurements of Fig. [2]). The same nomenclature applies to 
functions. In the case of Y-junctions, we denote them using the 
index m, I to specify the beam combination to which they are re- 
lated, with m corresponding to the actual beam that enters the 
Y-junction. The index k designates its two outputs. For exam- 
ple, Y[ 2 and Yj 2 are the two outputs of the Y-junction that splits 
beam 1 in signals that will interfere with beam 2. Similarly, Yl, 
and Y?j are the outputs of the Y-junction that splits beam 2 into 
signals that will combine with beam 1 . We use the same notation 
for the couplers, e.g. C\ A and C\ A are the outputs of the coupler 
corresponding to the combination of beam 2 and 4. 

With this notation, the intensity recorded at the outputs of 
the combination of beams m, I can be written : 



where N m is the number of photons in the m beam and t k ml the 
total transmission of the k output for the ml beam pair. V k , is the 
instrumental contrast, <p k , is the instrumental phase introduced 
by the IO beam combiner between the two interfering beams. 
ip p , is the residual atmospheric phase due to piston effects. V b * 




Fig. 3. Principle of the phase-shifting function: a variation of the 
optical path is induced by a differential change of the waveguide 
effective index due to a change in their width. The concatenation 
of carefully optimized portions of waveguides with controlled 
effective index allows to flatten the wavelength response. 



3. Laboratory set up 

In this section, we present the aims of the experiments, our 
testbed as well as our operating mode for the data acquisition 
and processing. 



3.1. Goals of the characterization 

With such beam combiners, all the information about the coher- 
ence of the object is included in the way the 4 pixels are related 
to each other, including the instrumental contribution. This con- 
tribution has therefore to be known, i.e. fully calibrated. 

The relationship between the measured fluxes on the pix- 
els and the visibility amplitudes and phases of the object can 
be expressed with a matrix representing the behavior of the in- 
strument. With an unresolved internal source (i.e. V ob J = 1 and 

v ml 

(p°J = 0) and without piston (i.e. centered at zero OPD), Eq. ([D 
becomes : 



(2) 



with X k , — Jt k ,t k cos(ui k ,), a coefficient different from 1, that 

ml y ml Im v ~ ml" ' 

corresponds to the level at which the beam combiner conserves 
the coherence and that depends on iff., the internal IO phase 
specific to the output ml k . If one isolates a combination cell [ml ], 
the relation between the output intensity and input number of 
photons (Eq. (0) can be written as : 



i \ 

ml 

f 
ml > 



( t 1 , t) V 1 ,X> ,(</>' ,) , v 

ml un ml ml ' ml' 

t 2 t 2 V 2 Y 2 (in 2 \ 

'ml l bn y ml A mlW m l> 

f , f, V i ,X i ,(w i ,) 

n^l Im r^l v ' ml* 



y/N m Ni ) 



(3) 



is the object visibility and ip ob J is its phase. 



ml 



ml 1 and ml 2 correspond to two outputs of the same coupler (the 
same is valid for outputs ml 3 and ml 4 ). Therefore, ideally, be- 
cause energy is conserved at the output of a coupler the follow- 
ing relations should apply: <p^ = tp) nl + n and ip ml = <p 3 , + n. 
Similarly, the beam combiner is ideally designed to introduce a 
phase quadrature between the outputs therefore: (p\ — <p x m[ + jt/2 
and w 4 , = w 2 , + n/2. 

ml 'ml 1 

The overall behavior of the beam combiner can be general- 
ized in a matrix. When considering all combinations, the matrix 
should then be constituted of similar blocks of a [4 x 3] ma- 
trix with zero elsewhere. In reality, crossing terms appear both 
as incoherent and coherent contributions, and the actual outgo- 
ing intensities should be described using a general matrix of 
24 x 10 terms : 
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i A 



'12 '21 



V 1 X 1 



V 4 X 4 

^34^34 



#1 



N 4 



(4) 



This matrix, called the V2PM (visibility to pixel matrix), 
in accordance to previous w ork on multiaxial interferometers 
dTatulli & LeBouquin 2006), completely characterizes the in- 
strumental behavior of the beam combiner e.g. the transmission, 
the visibility, the phase relations and the parasite flux. The ul- 
timate goal of such a study w ould be to precisely estimate and 
calibrate it (lLacour et al . 2008). However, it is out of the scope of 
such a paper to present and discuss a full characterization of the 
global beam combiner matrix. We prefer to focus on character- 
izing the individual tricoupler functions, Y-junctions, couplers 
and phase-shifting devices described by Eq. © as well as the 
global routing of the incoherent flux (so-called crosstalk) inside 
the beam combiner. As it will be seen later, the crosstalk terms 
(related to the crossing terms in Eq. (|4]l) are sufficiently small, 
which justifies this approach. 

To reach this goal, we set up a laboratory testbed and tested 
the IO beam combiners through photometric and interferometric 
measurements to calibrate this instrumental matrix. The follow- 
ing quantities, as well as their dependence on wavelength, have 
been measured: 



ml 



the so called "normalized kappa matrix" k 
instrumental contrast V*,; 

instrumental phase shift between outputs, ^Lr^t 2 , with (fr, 
being the individual phase of the output signal ml k ; 



4. instrumental closure phase O 



mlj 



<f k , 

'ml 



+ (fl*. + US k . 



jm 



3.2. Testbed description 

We designed a dedicated interferometr ic testbed c apable of sim- 
ulating an 8 telescope interferometer ( Jocou 2007). Figure|4]de- 
scribes individual functions of the setup. The bench includes var- 
ious items : 

(a) an object simulator that can reproduce a single star or a 
binary star with an adjustable flux ratio 

(b) up to 8 optical devices simulating telescopes and cou- 
pling the light into single-mode polarization-maintaining 
fibers 

(c) optical path compensation and modulation devices (delay 
lines of a few mm long) 

(d) an IO beam combiner 

(e) a spectrometer 

(f) a Wollaston prism to split the linear polarizations 

(g) an infrared detector. 

All laboratory tests were carried out in the H band with light 
sources of different coherence lengths. The object simulator can 
reproduce a single star as well as a binary star. In the latter case, 
its design is based on an optical setup that mimics a Michelson 
interferometer, but with a tilted mirror in one of its arm and with 
an unbalanced pathlength between the two arms. It produces 
two non-coherent luminous spots simulating a binary star, whose 
separation can be adjusted by tilting the mirrors. This setup will 
be used to characterize the dynamics of the testbed. The image is 
placed at the focal plane of a F/5 collimator to produce a 100 mm 



Telescope array 
(4 telescopes) 




Control/software 



Fig. 4. Schematic view of the laboratory testbed simulating the 
VLTI (see text for details). 
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24 outputs 



Fig. 5. A detector image of the 24 outputs, obtained when using 
the spectrograph and a Wollaston prism (splitting the two linear 
polarizations PI and P2). The patterns are due to various non 
zero OPD for the different beam combinations. 



diameter collimated beam. The wavefront is sampled by up to 8 
telescopes that can be set to reproduce a replica of the VLTI en- 
trance pupil. These telescopes are made up of an F/5 gradium 
lens feeding a polarization maintaining fiber and are placed in 
the collimated beam. Each telescope is mounted on an individ- 
ual module with tip-tilt adjustments and a motorized translation 
stage (delay lines). Shutters are used to block the light of each 
telescope. The fibers, that have been equalized to limit differ- 
ential effects (AL=lmm), are gathered in a V-groove chip that 
feeds the IO beam combiner. A diffracting grating provides 15 
spectral channels through the H band, while the Wollaston prism 
splits the linear polarization states to improve the instrumental 
transfer function. These two last elements, that can be placed 
and removed easily depending on the need, are located in an 
afocal mount with a magnification of 1 . The detector is a near- 
infrared InGaAs PICNIC chip, with 40/mi-large pixels. In the 
case of wide band measurements, each waveguide output is im- 
aged on one single pixel. Figure [5] is a detector image obtained 
when both the spectrograph and the Wollaston prism are used. It 
shows the 24 beam combiner outputs spectrally dispersed along 
the vertical direction. PI and P2 correspond to the two linear 
polarization states. 
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Sh = shutter; 0=open; X=closed 

Table 1. Experimental protocol including photometric and in- 
terferometric measurements. 



3.3. Data acquisition and processing 

3.3.1. Protocol 

The sequence of acquisitions performed to characterize the beam 
combiners consists of 6 steps that can be done with or without 
spectral dispersion (see Table [TJ. Step 1 is a background mea- 
surement with all shutters closed to prevent any light from prop- 
agating through the instrument. The 4 consecutive steps are mea- 
surements with only one beam at the time that give access to 
the flux splitting ratios in the couplers and tricouplers. Finally, 
Step 6 is the interferometric combination of all input beams. All 
measurements are repeated 1024 times. 

The first 5 steps are used to validate the design in terms of 
photometry (light routing, transmission and splitting ratios, un- 
desired flux, together with their wavelength dependence). Step 6 
leads to the determination of the value, stability and chromatic - 
ity of instrumental contrasts and closure phases as well as of 
the phase relations between phase-shifted outputs, supposedly in 
quadrature. To get a complete and independent laboratory char- 
acterization of each output of the tested beam combiners, the in- 
terferometric measurements presented in this paper are obtained 
with OPD modulation and polarization splitting on a point-like 
source. 

3.3.2. Data reduction 

The data processing derives four quantities : the kappa matrix 
(i.e. the photometric contribution of each beam to the interfero- 
gram); the instrumental contrast; the phase-shift induced by the 
devices and the closure phases. 

The kappa matrix is extracted from each individual set of 
data where only one input is illuminated (TableQ] Steps 2-5). 

The instrumental contrast is computed using a classical visi- 
bility estimator on the interferograms to evaluate the coherence. 
This consists of estimating the envelope amplitude and calibrat- 
ing for the photometric inbalance between the interfering beams. 
Since our experimental data are obtained with a SNR of ~100, 
the use of such a simple estimator is appropriate. 

The phase-shift between two outputs is computed as the 
phase of the complex product of the Fourier spectra correspond- 
ing to the two signals. For two phase-shifted outputs related to 
the combination cell [ml] (written ml 1 and mP, i.e. outputs k — 1 
and k — 3 respectively), it is calculated as follows: 



Ttnl Tn 



= atan 



( y[F m p(v m ii)F* m p(v m0 )} 



^[F,„p (v m / i i)FJ |P (v mW )] t 



where 3,% stand for the imaginary and real parts respectively. 
FK(v m i t k) (here with k = 1, 3) is the intensity of the Fourier spec- 
trum of the signal ml k , taken at the maximum value (to which 
Vmijt corresponds). 



The closure phase is measured using a triplet of telescopes 
i.e. from three pairwise combinations. It is calculated as the 
phase of the bispectrum, that is the complex product of the cor- 
responding three Fourier spectra. Since each beam combination 
produces 4 phase-shifted output signals, there are 4 closure- 
phase signals per telescope triangle. For the telescope triangle 
[mln], the closure-phase derived from the output k can be writ- 
ten as: 



mlj 



atan 



K [F^v^Ff/vt^F^iv^ J 



The closure phase is computed with the constraint that the 
frequencies respect the closure relation (e.g. v m/ / t = v m u + v//,jt)- 
Definitions of F*,(v m yt), F*(vy tj and F* .(v m ^) are identical to 
the phase-shift case. 



The methodology for data reduction with spectral dispersion 
is identical. From all illuminated pixels of the detector (in a case 
such as in Fig. |5), we measure interferograms from which we 
derive the chromatic behavior of instrumental quantities. 



4. Results 

In this section, we present the results of the laboratory experi- 
ments obtained on a point-like source, in terms of flux through- 
put and routing, instrumental contrasts, phase-shifts, and closure 
phases. The results correspond to two different 10 chips manu- 
factured in the same wafer (called Chipl and Chip2 in all tables) 
in broad band, as well as with spectral dispersion for the second 
chip only. 

For the latter experiment, we report in all tables the aver- 
age values as well as the amplitude of the variation over the 
wavelength range (i.e. \X„ mx ^ - X m j„^\, for the instrumental 
quantity X). We refer to the latter as chromaticity in the text. 
Detailed studies that relate the performance to the IO design 
and simulations will be given in a following paper (Labeye et 
al. 2009, in prep.). All results are commented on in Sect. The 
notations used in the tables and figures are the same as defined 
in Sect. |2] 



4. 1 . Photometric measurements 

Transmission: When injecting 100 photons in one input, the 
transmission is the total number of photons detected at all out- 
puts. The overall transmission budget includes the coupling ef- 
ficiency from the telescope point spread function to the fiber, as 
well as the propagation losses inside the fibers and the IO chip 
transmission. The latter quantity is determined in broad band by 
using a fiber at each input and at each output, and the measured 
flux is normalized by a fiber-to-fiber transmission. The measure- 
ment gives about 65% in the H band for the transmission of the 
IO chip itself. The telescope-to-fiber coupling is ideally as 80% 
using a perfect circular pupil without central obscuration. The 
silicate fiber transmission in the H band is excellent (< 5dB/km), 
consequently, for the 2m fiber lengths that we are using, the cor- 
responding transmission is * 99%. The total throughput of the 
'fibers+combiner' is therefore * 64%. 

Flux routing and individual IO functions: For the following 
paragraphs, we use Step 1 to Step 5 (Table [TJ. Figure [6j left, 
gives the broad band photometric coefficients for the 24 beam 
combiner outputs obtained when shutters prevent three telescope 
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Fig. 6. Left: kappa matrix photometric coefficients obtained when the light is injected in one input at a time (F1..4) for Chipl. The 
stars and squares show results obtained on two different days of experiments (called Day 1 and Day 3). Pixels from 1 to 24 are the 
24 outputs of the beam combiner, identified as 12 1 to 34 4 in Table [2] Right: variation of the kappa matrix photometric coefficients 
with wavelength over the H band range, measured on Chip2. The four panels correspond to the light injection in one input at the 
time (e.g. F[ for the injection in input 1). Inside each of them, among the 12 outputs illuminated, only 3 are plotted (in full, dashed, 
and dotted lines), corresponding to each combination cell. 



Table 2. Kappa matrix photometric coefficients obtained when the light is injected only in the 4th input (F4). The first two lines 
correspond to the wide band experiments while the third gives the average value over wavelength, as well as the peak-to-valley 
amplitude over the wavelength range (chromaticity). Values are divided by 10 4 . Bold numbers, preceded by a star, indicate the 
illuminated outputs. 
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706±9 


1069+9 


681+10 


Chip 2 


0.5±3 


0.5±3 0.8± 3 0.5±3 


2±3 


3±3 1±3 


1± 3 


2±3 


3± 3 


801±4 


932+4 


751+3 


avg / A/t 


1/2 


2/6 1/4 2/6 


3/9 


4/6 4/5 


3/5 


5/5 


6/11 


621/344 


926/357 


624/237 


Output 


*14 4 


*24 l *24 2 


*24 j 


*24 4 


23 3 


23 4 


*34' 


*34 2 


*34 j 


*34 4 


Chip 1 


1078±10 


800±10 1055±10 


811+ 10 


1037± 9 


4+10 


10±9 


710±9 


865±9 


675+ 9 


834+9 


Chip 2 


944±4 


784±4 930±4 


691 ±4 


990±4 


4+3 


3 ±3 


693±3 


896±4 


691+4 


874+4 


avg / AA 


953/184 


794/270 1032/150 


820/171 


1055/124 


8/6 


6/8 


585/194 


865/412 


766/360 


913/573 



beams from propagating through the chip. These coefficients are 
defined, for each pixel, by the ratio between the flux detected on 
one pixel and the sum of the flux on all outputs. As expected in 
such a case, 12 out of 24 pixels are illuminated. Two measure- 
ments taken on 2 different days are compared (stars and squared 
symbols) showing very small time variability. Table [2] gives an 
example of the averaged photometric coefficients for all 24 out- 
puts, obtained when the light is injected in the 4th input. 

With the same experiments using spectral dispersion, one 
can derive their dependence on wavelength. Figure[6] right, gives 
their variation with wavelength over the spectral range. Table [2] 
also gives the chromaticity of photometric coefficients obtained 



when injecting in the 4th input. 

To be more general, Table [3] gives measurements in all cases, 
i.e. when the light is injected in all 4 inputs, one at a time. We 
only provide the values corresponding to the outputs presenting 
the minimum and maximum chromaticity as well as the average 
chromaticity over the 12 signals. Because there is an important 
spread across the photometric coefficients values (see Table 
the chromaticity is given with respect to the coefficient value ob- 
tained when averaging over the spectral band (i.e. divided by this 
value). 

From these coefficients, in both wide band and spectrally dis- 
persed experiments, we determine the splitting ratio of the dif- 
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over 1024 measurements and is smaller than the symbol sizes. The theoretical values of 33% for the tricoupler, and 50% for the 
Y-junctions, are given by the horizontal dashed line. Values for all functions are given in Table [4] and [6] 



Table 4. Tricoupler splitting ratio measured in wide band. The dispersion (rms) over 1024 measurements is 0. 1 %. The two first lines 
correspond to the wide band experiments while the third gives the average value and the variation amplitude over the wavelength 
range. 



T 


T, 


T 2 


T 3 


T 4 




Chip 1 
Chip 2 

avg /A/l 


33.7 -33.3 - 32.8 
31.9-35.3 - 32.7 
32.2/10.5 - 37.2/2.9 -30.5/11.9 


33.2- 32.6- 34.1 
33.8- 34.0- 32.1 
32.9/18.6- 31.4/6.2- 35.6/13.2 


34.5 - 35.0 - 30.3 
32.4- 36.4-31.1 
32.1/2.6- 37.0/17.9 - 30.8/19.0 


34.2 - 35.9 - 

34.3 - 34.0 - 
31.3/9.9-37.1/5.0 


29.8 
31.6 

-31.4/14.2 



Table 3. Minimum, maximum and average of the photometric 
coefficients chromaticity among the 12 illuminated outputs from 
each injection. 



Injection in input # 


Average 


Minimum 


Maximum 


1 


30% 


14% 


48% 


2 


51% 


30% 


81% 


3 


58% 


27% 


80% 


4 


35% 


12% 


63% 



Table 5. Minimum, maximum and average chromaticity of tri- 
coupler splitting ratio (among the three outputs of each tricou- 
pler). Ti is the tricoupler corresponding to input 1. 



Tricoupler 


Average 


Minimum 


Maximum 


Ti 


8% 


3% 


12% 


T 2 


13% 


6% 


19% 


T 3 


13% 


3% 


19% 


T 4 


9% 


5% 


14% 



ferent optical functions (tricouplers, Y-junctions and couplers), 
under the assumption that all functions are ideal (z'.e no photon 
loss; Yti x i - 100%, with x-, a splitting ratio coefficient). For the 
sake of clarity, for each Y-junction and coupler, only one value 
out of the two splitting ratio coefficients is given in the tables, 
since the second output is obviously its complementary to 100%. 



Tricouplers: Table |4] gives such values for the 4 tricoupler 
showing flux splitting ratios close to 33% for both chips, in wide 
band. The best flux separation is 33.7/33.4/32.8, for the three 
outputs, with a rms over 1024 measurements of 0.1%. These 
values are similar a few days later with a variation from 0.1% 
to 1.5%. 

The results are comparable with the spectral dispersion. In this 
case, the closest splitting ratio from 33% is 32.9/31.4/35.6 + 
0.1%. Figure|7] left, presents the variation of the tricoupler split- 
ting ratio with wavelength. Table [5] gives the minimum, maxi- 
mum and average chromaticity of the tricoupler splitting ratio, 
among the 3 outputs of each of the four tricouplers of Chip2. 



Y-junctions: Table [6] gives the splitting ratio for the 12 Y- 
junctions. The values, measured in broad band, are close to 50%, 
with a 0.1%-dispersion over 1024 points. Variations are small 
from one day to another (2.9% maximum). 
With spectral dispersion, the splitting ratio are similar to the 
broad band measurements. Figure|7] middle and right, gives their 
wavelength-dependence. For clarity, only extreme behaviors are 
shown in the figure, with the smallest (middle plot) and greatest 
(right plot) variations over the H band. Out of the 12 Y-junctions, 
10 show a maximum variation inferior to 9% over the spanned 
range of the H band, while 2 show a strong variation of about 
26 and 28%. These two Y-junctions are the closest to the inputs 
(combination [23]). Over all the Y-junctions, the average chro- 
maticity is 9.2%. 
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Table 6. Y-junction splitting ratio. The dispersion (rms) over 1024 measurements is 0.1%. The two first lines correspond to the 
wide band experiments while the third gives the average value and the variation amplitude over the wavelength range. 



Y 


Y 1 

1 12 


Y 1 

1 12 


Y 1 

1 13 


Y 2 

1 13 


y' 

1 14 


Y 1 

1 14 


Y 1 

1 14 


x 24 


Y 1 

1 23 


Y 2 

1 23 


Y 1 

x 34 


Y 2 

x 34 


Chip 1 


52.5 


47.4 


47.0 


44.4 


50.2 


50.2 


52.6 


50.0 


51.5 


52.6 


54.7 


51.0 


Chip 2 


48.7 


45.7 


52.2 


49.8 


51.2 


50.5 


54.1 


50.4 


53.0 


53.2 


52.2 


50.3 


avg / AA 


52.1/4.0 


46.6/6.0 


52.0/5.7 


48.8/8.0 


49.7/7.5 


49.3/7.4 


51.6/3.7 


49.3/4.2 


49.5/28.0 


50.1/25.9 


49.5/4.5 


46.5/6.2 
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Fig. 8. Variation of coupler splitting ratio with wavelength. Only one of the flux ratio coefficients is plotted, and extreme behaviors 
are given : left, when the light is injected in the first input; right, for an injection in the second input. The theoretical value of 50% 
is given by the horizontal dashed line. 



Table 7. Coupler splitting ratio. The dispersion (rms) over 1024 measurements is 0.2%. The two first lines correspond to the wide 
band experiments while the third gives the average value and the variation amplitude over the wavelength range. 



c 


C 1 

12 


C 2 


C 1 

^13 




C 1 

14 


c- 

14 


C 1 


c- 


C 1 

1 -23 




C 1 

^34 


C 2 
^34 


Chip 1 


54.1 


59.6 


50.5 


58.4 


59.7 


61.5 


57.5 


58.5 


48.3 


50.4 


53.9 


55.2 


Chip 2 


52.0 


58.9 


55.2 


55.7 


58.3 


56.3 


56.3 


51.2 


52.6 


55.0 


55.0 


55.4 


avg / AA 


53.9/5.4 


53.9/4.1 


55.6/4.5 


57.2/5.7 


56.0/6.1 


60.0/5.8 


56.0/3.1 


56.5/5.4 


54.8/7.0 


47.6/5.0 


47.2/5.0 


54.1/7.2 



Couplers: Table [7] gives the flux splitting ratio given by the 
12 couplers, in broad band, showing asymmetric splitting, up to 
61.5%/38.5%, similarly for the two beam combiners. 
With spectral dispersion, the obtained values are similar. The 
variation of the 12 coupler splitting ratio with wavelength is 
given in Fig.[8](with only extreme behaviors). All couplers show 
low chromaticity, with an average maximum variation of about 
5%, with minimum and maximum values of about 3 and 7%. 

Finally, although, in theory, all functions are identical, we 
notice some slight variations from one another. In broad band, 
for Chipl, the differences are up to 1.7%, 2.7% and 4%, for 
the tricouplers, Y-junctions and couplers, respectively, while for 
Chip2, they are of 0.9%, 2.3%, and 2.6% for the same functions. 

Cross-talk: On the pixels where one should not detect any 
flux, the measured intensity gives the amount of undesired flux, 
i.e. cross-talk flux. It can be due to direct propagation into the 
substrate and to light leak at the X-junctions level, where waveg- 
uides are crossing. We estimate the cross-talk flux to be less than 
1.2%+0.4% of the total flux. Table [8] gives the measured cross 
talk flux when the light is injected in one input at the time. The 



impact of such an amount of cross-talk on the photometric quan- 
tities is within their error bars. Also included in Table [8] is the 
average cross-talk over wavelength with its chromaticity. 

4.2. Interferometric measurements 

Instrumental contrasts : Figure [9] left, shows an example of 
the interferograms obtained with Chip 1 in broad band. From 
these interferograms, we derive the instrumental contrasts af- 
ter calibrating for the photometric inbalance between interfering 
beams. Table [9] gives the measured values, showing high con- 
trast values. The minimum and maximum values are respectively 
of 95%+l% and 98%+l% for the first beam combiner, and of 
82%+6% and 94%+ 1% for the second one. 
These results show a maximum variation of about 5% over a 
day timescale and 10% from one day to another. The measured 
contrasts on the four phase shifted outputs are only slightly dif- 
ferent (on average 2%, up to 9%). The non-perfect linearity of 
our detector can lead to a bias of up to 5% in visibility. In the 
case of our experiments, this effect could not be reproduced and 
calibrated. Therefore, although the statistical errors can be very 
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Fig. 9. Left: an interferogram obtained with temporal OPD modulation with wide-band measurements. Right: the variation of 
instrumental contrasts with wavelength, for all four outputs of each beam pair. 100% contrasts are given by the horizontal dashed 
line. 



Table 8. Cross-talk flux (in %) determined when the light is injected in the 4 outputs successively. The two first lines report results 
obtained in broad band while the third line corresponds to the spectral dispersion case. Error bars are estimated from the dispersion 
(rms). 



Day 


Injection in 1 


Injection in 2 


Injection in 3 


Injection in 4 


Chip 1 


0.7±0.2 


0.6±0.2 


1.2+0.4 


1.1+0.3 


Chip 2 


0.4±0.2 


0.3+0.2 


0.3+0.2 


0.2+0.2 


avg / AA 


0.3/0.03 


0.4/0.03 


0.9/1.5 


0.5/0.6 



Table 9. Instrumental contrasts (in %) obtained with wide band experiments (first two lines) and with spectral dispersion (third 
line) with a statistical error of 1%. 



Output 


12' 


12 1 


12 J 


12 4 


23 1 


23 2 


13 1 


13 i 


13 J 


13 4 


14' 


14 a 


Chip 1 
Chip 2 

avg / AA 


95+1 
90+2 
99/2.7 


97+1 
90+2 
100/2.1 


95+1 
92+2 
99/2.6 


98+1 
88+3 
100/1.8 


96+1 
85+3 
98/3.3 


97+1 
93+2 
99/3.8 


94+1 
83+5 
97/1.0 


96+1 
88+4 
97/1.3 


94+1 
85+4 
97/0.8 


95+1 
88+3 
97/1.4 


96+1 
88+4 
97/1.3 


98+1 
94+1 
97/1.6 


Output 


14 J 


14 4 


24 1 


24 z 


24 J 


24 4 


23 J 


23 4 


34' 


34^ 


34 J 


34 4 


Chip 1 
Chip 2 
avg / AA 


97+1 
88+4 
97/1.2 


97+1 
93+1 
97/1.5 


96+1 
84+4 
100/2.5 


97+1 
93+2 
98/2.4 


97+1 
88+3 
99/2.2 


97+1 
91+2 
97/2.4 


97+1 
86+4 
98/1.1 


97+1 
92+2 
99/2.6 


96+1 
82+6 
97/3.0 


97+1 
88+3 
97/4.1 


97+1 
83+8 
97/3.5 


96+1 
84+3 
97/4.6 



small (~1%), an additional bias of 5% affects the contrast val- 
ues. 

With spectral dispersion, the measured contrasts are very high 
(up to 100%) showing very small variations with wavelength 
(see Figure [9] right). The maximum and minimum chromatic- 
ities, among all 24 outputs, are of 4.6% and 0.8% respectively, 
with an average of 2.3%. 

Phase shifts: Figure[10]shows 4 phase-shifted interferograms 
in each panel, that correspond to the intensity of the 4 phase- 
shifted outputs for the interferometric couple [34]. These 4 out- 
puts are in different phase states, as it can be seen, and in these 
specific examples, in the left panel, the phase-shift is close to the 
expected value of quadrature. 



Table 12. Examples of phase-shifts (in degrees) measured in the 
two polarization states (PI, P2), for the 6 beam pairs. Statistical 
errors are of 1°. 



Beam Pair 


[12] 


[13] 


[14] 


[24] 


[23] 


[34] 


PI 


62 


87 


34 


67 


77 


88 


P2 


55 


80 


37 


62 


70 


77 



Table [10] gives the values of the phase-shifts obtained for the 
outputs designed to be in quadrature. For the first chip, on 5 of 
the 6 interferometric combinations, the phase shifts are close to 
the quadrature (from 78°±1° to 82°±1°). For the sixth phase- 
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Fig. 10. The phase-shifted interferograms recorded for the 4 outputs (full, dashed, dotted, dot-dashed lines) of the [34] beam pair 
are close to quadrature (left). On the contrary, the 4 outputs of the central combination [14] produce interferograms only slightly 
phase-shifted (-26°) (right). 
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Fig. 11. Left: wavelength dependence of the phase-shifts for the outputs designed to be in quadrature, for all 6 pairwise combinations 
([ml]). Error bars on each spectral point are of 1°. The theoretical value of 90° is given by the horizontal dashed line. Right: variation 
of the closure phase over the H band, for the triangle [134]. The four symbols correspond to the four outputs in quadrature. 



Table 10. Phase-shifts obtained with the wide band experiments (two first lines) with errors of 1 .0° (the dispersion (rms) over half 
a day); The third line gives the average values over the wavelength range. 



Phase-Shifts 


#12 O 


#13 O 


Ol4 O 


* 24 O 


4>23 O 


O34 O 


Chip 1 


79.0 


81.5 


26.1 


79.4 


81.6 


77.9 


Chip 2 


62.5 


87.3 


32.9 


67.5 


77.5 


87.7 


avg / AA 


71.6/4.1 


75.7/3.1 


29.5/3.2 


72.8/5.8 


77.1/13.0 


69.6/4.3 



Table 11. Closure-phase measurements for one independent triangle, for Chip 2 with spectral dispersion, as well as amplitude of 
the variation over the H band range. Statistical errors are of 2.5°. 



Triangle [134] 




«f« O ®B4 O 


®k O 


avg / A/i 


114.6/22.8 


-68.2/18.1 -1.9/22.2 


179.0/30.1 
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Table 13. Contrasts and phase-shifts obtained when the split- 
ting of polarizations is done before or after the combination, and 
when no splitting is achieved. Errors are statistical. 



Pixel 


12 1 


12 1 


12 J 


12 4 


< V > before 


83.7±0.3 


83.6±0.3 


81.3+0.3 


81.7+0.3 


< V > after 


79.6±0.2 


79.4±0.2 


78.2+0.2 


77.2+0.2 


< V > 


71.3+0.5 


71.2+0.5 


66.7+0.5 


66.9+0.5 



Polar splitting 


before 


after 





Phase-shifts 


78.98+0.03 


79.80+0.08 


76.12+ 0.07 



shifting function, corresponding to the central [14] combination, 
the measurement gives about 26°+ 1 ° . The right panel of Fig. [TOl 
corresponds to such a combination. The second chip gives dif- 
ferent results, with phase shifts spanning a larger range of values 
(e.g 62°±1°; 88°±1°). The function corresponding to the [14] 
combination still shows a much smaller value of about 33° + 1 .0°. 

Over a timescale of half a day, measured phase shifts show 
variations of 1° at most, and from one day to another, a maxi- 
mum variation of 3°. 

With spectral dispersion, the obtained values for the phase 
shifts are similar. Figure QT] (left) and Table [TOl show the wave- 
length dependence of the measured phase-shifts for the outputs 
expected to be in quadrature. The maximum chromaticity is of 
13° (for the [23] combination, as for the Y-junctions) and the 
minimum is of 3.1°. Over all outputs, the average chromaticity 
is 5.6°. 

Closure phases: Table QT| gives a set of independent instru- 
mental closure phases measured for the 4 phase shifted outputs 
of one telescope triangle, with spectral dispersion. The measured 
closure phases have non-zero values which means that the beam 
combiner functions themselves contribute to the phase budget. 
These terms can result from additional OPD originating in small 
length differences between waveguides or from the delay in read- 
ing the detector pixels. The phase relation between the vari- 
ous beam combiner outputs can be found again in the closure- 
phase values. In fact, the closure-phases differ by about 180° for 
couplers outputs in phase opposition and similarly, the closure 
phases measured at outputs theoretically in quadrature are dif- 
ferent to the corresponding phase sum between the telescopes 
(i.e. <p m i + (pij- <Pmj)- 

Figure Qj] (right) and Table QT| give the variation of the in- 
strumental closure phases with wavelength, in the case of one 
triangle of telescopes ([134]). The minimum variation over the 
wavelength range is of 18° while the maximum is of 30°. 

Polarization: The two linear polarization directions are per- 
pendicular to each other following the symmetry axis of the chip 
itself (i.e. within the beam combiner plane (horizontal) and par- 
allel to the light wavefront (vertical)). Such orientations are de- 
fined at the time of manufacturing the chip, and were confirmed 
by laboratory measurements. 

In our study, disparities between measurements obtained on 
both linear polarization states were noticed. Contrasts and phase- 
shifts can be up to, respectively, 13% and 10° different, meaning 
that the two polarizations propagate differently. Also, for one 
polarization state, instrumental contrasts are higher than for the 
second one for all beam combiner outputs except for the ones 
related to the central phase shifting function ([14] combination). 
Similarly, phase shifts are closer to the quadrature for one po- 
larization than the other for all outputs except the central ones 



(Table IT2b . Therefore, with respect to polarization, one phase 
shifting function has a different behavior than the 5 others. 

In order to identify possible problems with instrumental po- 
larization and to determine the best instrumental set-up to reduce 
the visibility loss due to birefringent fibers, a detailed study on a 
day timescale was done with three optical set-ups: the first one 
had polarizers before injection into the fibers and beam com- 
bination (i.e. before the telescope mounts); the second one had 
a Wollaston prism before imaging on the detector, after beam 
combination; and the third one had no polarization splitting. 

We found that the instrumental contrasts can drop by 10 to 
15% when no polarization splitting is done (Table [T3l. Slightly 
better contrasts (~3%) were found when using polarizers be- 
fore beam combination rather than after with a Wollaston prism. 
When comparing stability for all configurations, no significant 
difference was found. 



5. Discussion 

The characterizations presented here shows that global proper- 
ties of the designed beam combiners are very satisfactory for a 
first prototype. In this section, we discuss how departures from 
the ideal case might affect the performance. 

Photometry: The transmission of the beam combiners di- 
rectly impacts the instrument sensitivity and therefore, con- 
strains the limiting magnitude. For these longest and most com- 
plex IO beam combiners tested today, a 65% transmission is sat- 
isfactory. An improved technology has allowed to reduce the 
losses with respect to previous beam combiners. For compar- 
ison, the IONIC3T/IOTA H band beam combiners that were 
made with the first technology and had only 3 Y-junctions and 
3 couplers, presented a transmission of 60%. The gain comes 
essentially from improvement in propagation losses associated 
with a reduced beam combiner length. 

We show evidence for small crosstalk photon leaks that lead 
to unwanted flux in the combining cells. These potentially affect 
the photometry estimation and might introduce a small coher- 
ent contribution not revealed in this study. However the mea- 
sured values for both chips are mostly inferior to 1 % of the total 
flux and have contributions smaller than the typical error bars 
of the measurements. The origin of this effect has been identi- 
fied as coming from imperfect fiber/IO coupling (in the exper- 
iment the fibers are not glued unlike in an actual instrument) 
and from flux coupling in the substrate that is partially guided. 
However, the latter contribution has been dramatically reduced 
thanks to new etching technology that allows each waveguide 
to be completely iso lated from the others and to not waste flux 
dLabeve et alj |2006). As a matter of fact, our experiments using 
the old technology showed flux crosstalk of up to 8%. 

The flux routing of the interfering beams, for each pair, needs 
to be as equal as possible, although this is not a strong require- 
ment. In fact, the instrumental visibility (and therefore the SNR) 
decreases as the beam fluxes are unbalanced. For these beam 
combiners, the tricouplers equally split the flux in three (within 
the error bars), and the flux splitting ratios of the Y-junctions are 
satisfactory. On the contrary, the couplers can be as unbalanced 
as 60%/40% leading to a maximum contrast loss of 2%. These 
splitting ratios were estimated with the assumption of ide al func- 
tions, which actually present a small loss of a few % {Cabeve 
2008). These unbalanced ratios are due to an error in the de- 
sign, and new prototypes are being made with improved cou- 
plers. Globally, the agreement of the simulated data with the ex- 
periments is very satisfactory. 
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Finally, the stability and repeatability of all observables are 
key elements to reach high dynamics. Photometric quantities are 
stable over a day timescale at the level of a percent, and are re- 
producible from one day to another with a maximum variation 
of 3%. We suspect that the absence of a glued interface between 
the fibers and IO beam combiners contributes to the small errors 
(except for the closure phases for which this effect is canceled 
out). However, for a first demonstration, the performances are 
good and the general flux routing is validated. The competitive 
transmission together with the capacity to use all photons for 
coherent combination lead to an overall high sensitivity for such 
4-beam phase-shifting beam combiners. 

Interferometry: The instrumental contrast directly impacts 
the sensitivity of the instrument and should therefore be as high 
as possible. Here the beam combiners produce instrumental con- 
trasts with values always above 80% in wide band, and reaching 
100% with spectral dispersion. In broad band operation, the ef- 
fects of differential dispersion due to unequal length of the fibers 
can lead to a contrast loss of up to 20% which explains the very 
high contrasts obtained with spectral dispersion. Statistically, the 
stability of the contrasts over a day timescale is good (from 1% 
to 4%) but the reproducibility of such measurements onto an- 
other day could not be validated due to the non-linearity of our 
detector. This affects the value with a bias of 5%, making our 
results only upper limits. 

The phase shifting functions were initially designed to sam- 
ple the coherence in four phase states in quadrature (the so called 
ABCD sampling). For the two tested beam combiners, the phase- 
shifting functions lead to 5 phase-shifts out of 6 in agreement 
with our expectations, while one of them (the central one corre- 
sponding to the [14] combination) is far from 90°. We suspect an 
inhomogeneity in the constitution of the silicon substrate at the 
position of this phase shifting function. Of all 6 functions, this 
one is the most distant from the center of the beam combiner 
and is located next to the edge of the chip. This effect could ex- 
plain the measurements obtained on both chips since they were 
located next to each other on the same wafer. This could also 
result from a relaxation of the stresses on the silica during the 
chip cut. An error in the design has been ruled out. This is the 
first time that these functions have been tested and while still not 
completely controled, these results are encouraging and show 
that the use of phase shifting functions in IO beam combiners is 
promising. In addition, phase shifts are stable on a timescale of 
a day, within 1° and vary up to 3° from one day to another. 

The incidence of this departure from phase quadrature on 
the final complex visibility SNR estimation cannot be quantified 
without a proper calculation. It should be seen as a reduction in 
the instrumental response. In the limiting case where the phase 
shift is 0, one cannot retrieve unambiguously the complex vis- 
ibility information. However since all but one phase shift have 
values close to quadrature, we believe this validates the concept. 

Finally, the measured closure phases are not equal to zero. 
Specific phase contributions of each function, that result from 
non-perfectly symmetric pathways, are unknown. When observ- 
ing a scientific target, calibrating with a point source (that is 
centro-symmetric and supposedly leads to a zero closure-phase) 
should allow us to remove the instrumental contribution. It 
should be expect ed that given the re markable stability of such 
beam combiners (Bergeret al. 2000) the instrumental response 
should be very well calibrated. 

Chromaticity : The importance of the chromaticity of the 
functions directly depends on the spectral resolution used in the 
instrument. 



On all photometric coefficients, the average chromaticity can 
go from about 12% to 80%. Individual functions show varieties 
of chromaticity, with variations for tricouplers and Y-junctions 
from 3% to 28% while couplers can be considered as achro- 
matic. The Y-junctions and couplers for the [23] combination 
are more chromatic than the others. More details will be given in 
a following paper (Labeye et al. 2009, in prep.). However photo- 
metric calibration in the presence of dispersion should solve this 
issue. With correction for photometric effects, the obtained con- 
trasts show very little chromaticity, with an average maximum 
variation of 4.6% over the range of H band. Phase shifts show 
chromatic variations of 5.6° on average, the maximum variation 
corresponding to the [23] combination as for individual func- 
tions. The closure phases show strong variation with wavelength, 
up to 30°, which results in the sum of all chromatic effects for 
the three telescopes in the closure relation. 

We can anticipate that in cases where spectral resolution is 
low the effective wavelength might be affected and a proper 
wavelength calibration should be considered together with a 
proper stellar calibrator choice. In addition, in order to limit 
potential biases, particular care should be taken to specify the 
alignment accuracy between the beam combiner, spectrograph 
and detector pixels. 

Polarization: Birefringence control is a critical part of a 
guided optics instrument. In our experiment, we used highly 
birefringent fibers which have a well defined polarization axis 
but, in turn, require specific care on how the polarization state is 
modified along the propagation. 

Our study first shows that all estimated quantities (photome- 
try, instrumental contrast, phase shift, closure phase) can be dif- 
ferent depending on the polarization state. We have shown in 
particular that the behavior of the central [14] combination is dif- 
ferent from the other phase-shifting functions. This confirms the 
suspected marginal behavior of this combination, maybe origi- 
nating from an inhomogeneity in the substrate for both chips. 
Besides, our work showed that it is necessary to split the polar- 
ization states before or after the beam combinations, or to ac- 
tively control the phase shift between the two linear states to 
avoid contrast loss. The use of birefringent fibers and waveg- 
uides forces the phase velocity of the two polarization states 
to be different and results in two shifted interferograms. When 
these two are superimposed, it leads to a single low contrast in- 
terferogram. In the case of a sensitive imaging instrument, split- 
ting after the beam combination is recommended since the use 
of polarizers before the combination would lose half of the use- 
ful photons. As far as stability is concerned, without splitting, 
the differential phase shifting of the polarization states inside the 
beam combiner would lead to varying results. The new proto- 
types are being designed with special care given to this problem. 

On-sky operating modes: By allowing all the complex coher- 
ent factors to be measured in one single detector frame readout, 
these beam combiners offer two observational modes, depend- 
ing on the stability of the fringes (i.e. on the atmospheric con- 
ditions or on the availability of a fringe tracker and its perfor- 
mance). If the fringes are stabilized to better than a fraction of 
the wavelength, a long coherent integration of the flux on each 
pixel is possible (i.e coherencing mode), highly increasing the 
SNR compared to temporal encoding. Otherwise, by varying the 
OPD, one will access 4 phase-shifted interferograms on which 
to estimate the interferometric observables. The latter mode is 
also well suited for laboratory measurements and calibration. 

Precision interferometry and data reduction: 
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We have shown that the described beam combiners present 
performance well suited for astronomical beam combination in 
a four telescope imaging interferometer. Our extensive labora- 
tory characterization shows that on-sky performance, in terms 
of precision, should be comparable to what has been achieved 
with IONIC- VINCI/VLTI, IONIC3/IOTA or FLUOR/IOTA- 
CHARA. However, an interesting number of astrophysical prob- 
lems will soon be more demanding (e.g. debris disk, hot Jupiter 
detection). In that case it is important to better characterize 
calibration issues and tackle all imperfections that could be 
introduced by the beam combiner. Such work is justified by 
the tremendous stability of the beam combiner that has been 
revealed by numerous industrial developments. Therefore, a 
proper calibration of the beam combiner should allow system- 
atics biases to be removed. What we propose is to use the so- 
called 'visibility to pixel matrix' calibration in order to carry out 
a global inversion of the matrix that links the measured intensi- 
ties with the properties of the scientific object. Such a method, 
that will be detailed in a coming paper allows all instrumental 
contributions to be extracted. We briefly discuss here the philos- 
ophy of this data reduction. By developing the cosine, Eq. (Q]i 
can be written in such a way that instrumental terms are sepa- 
rated from the object contribution (V°^, <p° /, or similarly the 
complex visibility V), leading to a system of linear equations 
that includes the instrumental matrix (V2PM, see Sect. 

I = V2PM * V 

The complex visibility of the scientific source will be ob- 
tained by inverting the system. A full characterization of the 
V2PM matrix should be first done with internal calibration pro- 
cedures. This method allows all instrumental effects to be in- 
cluded in a single matrix taking into account all crossing terms 
due to unideal behavior. 

6. Conclusion 

We have presented a laboratory characterization, in the H band, 
of an integrated optics beam combiner dedicated to the combi- 
nation of a four telescope interferometer. It uses a novel 'pair- 
wise static ABCD' combination scheme which optimizes the ex- 
traction of coherence information: visibility amplitudes, phases 
and closure phases. Our measurements show that, although com- 
plex, the flux routing inside the beam combiner is efficient and 
that the global throughput is competitive (65%). In particular 
the comparison with simulated performance of several building 
block functions is very satisfactory. The instrumental contrasts 
are high (>80%) and stable showing maximum variations of a 
few % within a day. The instrumental closure phases are differ- 
ent from zero but stable and probably caused by small internal 
optical path differences. The specificity of these beam combin- 
ers, which is to produce four phase shifted outputs to sample 
the coherence information (for each baseline) in one integra- 
tion, has been validated. Finally, the global chromaticity of the 
beam combiner has required specific optimizations that are non- 
standard with respect to standard telecommunication 10 func- 
tions. IO chips used in telecommunications are usually required 
to show a flat response on short bandpasses (AA «20-30nm) 
while here we have extended the couplers and tricouplers and 
phase shifter r esponse to more than * 150 nm around 1.6/im 
(lLabevd l2008). These results on a first prototype validate the 
feasibility of a 'pairwise static ABCD' combination scheme and 
its suitability for an interferometric imaging instrument. 

Most of technological building blocks are now defined. We 
are working on the definitive version that will be included inside 



Gravity and VSI. This will require a number of technological 
improvements and innovations. In particular, the phase-shifting 
function will be ameliorated to get closer to the nominal 90° 
phase shift. The throughput will be improved with an optimized 
routing that will reduce the global propagation losses. We are 
currently extending the demonstration, using the same silica on 
silicon technology, to the K band, as required by Gravity and 
VSI. It is expected that, while less transmissive in this band, the 
short propagation distances inside the combiner will lead to ac- 
ceptable global losses. Finally, we will explore how this combi- 
nation concept can be extended to a six-way beam combiner and 
fit VSI requirements. 
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